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Abstract 

Results are presented of a study of the experimental and theoretical accuracy one may 
achieve at HERA in measuring the strong coupling constant and the gluon distribution 
from scaling violations of F2 structure functions. 



Introduction 

Accurate measurements of F2 structure functions in deep inelastic scattering provide one of 
the cleanest tests of perturbative QCD. In the first few years of experimentation at HERA the 
available kinematic range was extended to low values of x ~ 10~^ and large ~ 5000 GeV^. 
In the next 8 years of operation HERA might deliver integrated luminosities of 0.5-1 fb^^. The 
extended kinematic coverage and increased luminosity allows for detailed measurements of the 
scaling violations in F2 and hence of the strong coupling constant and the gluon distribution 
xg{x) at low x. 

In this report I will summarise the results of studies done during the workshop 'Future Physics 
at Hera' |l|, |^, ^ to estimate the experimental and theoretical errors on as(M|) obtained from 
a QCD analysis of future F2 structure function data. Also shown is the precision one may reach 
in the determination of the gluon distribution. 



Comparison of NLO evolution codes 



In this section I present results of a study on how well the various implementations of pertur- 
bative QCD are numerically and conceptually under control |l|. 

The evolution equations for the parton distributions f{x, Q"^) in the proton are given by 



dfix,Q' 
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asiQ^)Poix) + al{Q^)Pi{x) + 0{al) ® /(x, Q 



(1) 
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where we write as{Q'^) = as((5^)/47r and where Pq and Pi are the leading order (LO) and NLO 
sphtting functions respectively. 

Two methods are widely used to solve Eq. (|1]). In the first approach ('x-space') the parton 
distributions are numerically evolved on a grid in x and Q^. This method is conceptually simple 
but the numerical accuracy depends on the number of gridpoints in x and which is limited 
by the amount of CPU time one can afford to spend in the computations. 

In the second approach ( W-space') the Mellin transform of Eq. (|I|) is taken so that the convo- 
lution integrals become simple products. The resulting ordinary differential equations can be 
solved analytically. The result is then transformed back to x-space. This method is mathe- 
matically more involved but accuracies of ~ 10~^ are readily achieved @]. 

Detailed comparisons were made of two x-space programs used by the ZEUS (B) and 
HI (PZ) p collaborations and two A^-space progams which we label (V) [0 and (R) respec- 
tively. With these four programs identical sets of input parton distributions (10~^ < x < 1) 
were evolved in NLO from the input scale = 4 GeV^ up to = 10^ GeV^. 

Fig. |1| shows a comparison of the results. One notices the very good agreement between the 
iV-space programs (V) and (R). The two x-space programs (B) and (PZ) agree to within 
0.05% with the A^-space programs over a wide kinematic range. The agreement is slighly worse 
(~ 1%) at very high x where the parton distributions vanish. We remark that the agreement 
between the A^-space and x-space programs can in principle be further improved by increasing 
the number of gridpoints in the latter. 
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Figure 1: The relative difference between the up- valence (m^), singlet (S) and gluon (g) densities 
as obtained from evolving identical input at = 4 GeV^ in NLO with the evolution programs 
i = (B, PZ, V) and (R). 
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Experimental errors on cXg 



To investigate the experimental error on as{M'^) QCD fits were performed to simulated HERA 
F2 datasets listed in Table . The data cover a kinematic range 1.5 x 10~^ < x < 0.7 and 
0.5 < Q2 < 5 X 10^ GeV^. 

Seven independent sources of systematic error were taken into account (see for details) giving 
a total systematic error of ~ 1-5% over most of the kinematic range which is a factor of 2-5 
better than presently achieved. Residual systematic effects were represented by a point to point 
uncorrelated systematic error of 1%. 

The following model was fitted to the simulated data: 

F.{p,s) = Fp^''{p) (l-E^^^r) (2) 

I 

where F^^^{p) is the QCD prediction for F2, A^^f** is the (relative) systematic error on dat- 
apoint {€) stemming from source (/) and Si are the systematic parameters. It is assumed that 
these parameters are uncorrelated and gaussian distributed with zero mean and unit variance. 

The parameters (p) in Eq. (H) represent a<((M|) and those parameters describing the parton 
distributions at the input scale Ql = A GeV^. The gluon distribution {xg), the quark sin- 
glet distribution (xS) and the difference of the up and down quark distributions {x/S.ud) were 
parametrised as 

xh{x,Ql) = AhX^'^{l-xf'^P{x) (3) 

with P{x) = 1 for xg and xA^d and P{x) = 1 + D^/x + Ex for xS. In the studies presented 
below two types of fit were considered: (i) leave the parameters (p) and si free in the fit and 
(ii) fix the systematic parameters s; to zero. In the latter fits the systematic errors (As; = 1) are 
propagated to the covariance matrix of the fitted parameters (p) using the technique described 
1^. Since we are only interested in the errors the data were replaced by the model so that 



m 



the fits immediately converged. 

The results for the total error on are given in Table . A fit of the proton datasets I and 
II with a cut of of 3 GeV^ yields Aa,(Mf ) = 0.006 (0.012) depending on whether the 
systematic parameters are fitted or fixed (fit 1 in Table ). The error on is much improved 
when pertubative QCD is assumed to be valid at lower values of and the cut is lowered 
from 3 to 1 GeV^ (fit 2). Doubling the luminosity of the high sample has no effect which 
illustrates the fact that the error on as is dominated by the experimental systematic errors 
(fit 4). A modest improvement in the error is obtained when lower energy proton data are 
included (fit 6). 



Dataset 


nucleon 


Ee (GeV) 


En (GeV) 


L (pb-i) 


QLn (GeV^) 


Qla. (GeV^) 


I 


proton 


27.6 


820 


10 


0.5 


100 


II 


proton 


27.6 


820 


500 


100 


50000 


III 


proton 


27.6 


400 


200 


100 


20000 


IV 


proton 


15.0 


820 


10 


0.5 


100 


V 


deuteron 


27.6 


410 


10 


0.5 


100 


VI 


deuteron 


27.6 


410 


50 


100 


20000 



Table 1: Summary of simulated data sets for this study. 
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Fit 


datasets 


QLt (GeV) 


L(set II) (pb-i) 


syst. fitted 


syst. fixed 


1 


I, II 


3 


500 


0.006 


0.012 


2 


I, II 


1 


500 


0.002 


0.006 


3 




8 


500 


0.007 


0.012 


4 


i'ii 


3 


1000 


0.006 


0.012 


5 




3 


10 


0.010 


0.015 


6 


I, II, III, IV 


3 


500 


0.005 
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I, II, V, VI 


3 


500 


0.004 


0.009 




Table 2: The 


error Aa;s(M| 


) for several fits 


to Hera data 


alone. 



In the fits 1-6 xAud{x,QQ) was kept fixed to the nominal input value. The error on this 
distribution was taken from the QCD analysis of ref. and contributes 0.004 to the error on 
as(M|) (not included in Table ). This error is eliminated if HERA deuteron data are included: 
the difference of and constrains xA^d which can thus be left free in the fit. It is seen from 
Table that the error on is reduced even though the number of fit parameters has increased 
(fit 7). 

To investigate if HERA can improve the error on as from fixed target data fits were 
performed including SLAC 0, BCDMS ^ and NMC H proton and deuteron F2 data. To 
remove higher twist effects a cut > 10 GeV^ was imposed. A fit to SLAC and BCDMS 
data alone reproduced the result of ref. [l^: A«s(M|) = 0.003. 



Including the fixed target data the error on is much less sensitive to the cuts imposed 
and the luminosity of the high HERA sample. Depending on the cuts and the available 
luminosity the fits yield Aas(M|) = 0.0015 - 0.0020 when all systematic parameters are left 
free and Aa,(M|) = 0.0025 - 0.0035 when they are kept fixed. 



Of course, contrary to an analysis of high x fixed target data alone, the QCD fits to HERA 
structure functions result in a joint determination of both as and xg. For instance a fit to the 
datasets I and II (fit 1) with all systematic parameters left free results in a determination of the 
gluon distribution with an accuracy of about 3% at x = 10"'^ and = 20 GeV^, as illustrated 
in Fig. 



Theoretical errors on a. 



The following possible sources of theoretical uncertainty were investigated 

• The effect arising from different representations of as- The scale dependence of the strong 
coupling constant reads 

5a, (Q2) 



(91nQ2 



(3oaiiQ')-f3iai{Q')+0{at) 



(4) 



This equation can easily be solved numerically given a, at some input scale Q^. The 
following approximate solution in terms of the QCD scale parameter A is widely used: 

1 /3ilnln(QVA^^ 



/3oln(QVA2) (31 In^Q^A^) 
4 



+ 0(ln-^(Q7A^)) 



(5) 




Figure 2: The gluon distribution versus x for different values of obtained from a fit to 
simulated high energy HERA proton data. The bands show the total error. 

Either using Eq. (^) or Eq. causes a shift in q;s(M|) of 0.001 or less. 

• The offsets originating from the different prescriptions of the NLO evolution. In the A^- 
space approach the evolution equations and their analytic solutions are usually expanded 
as a power series in and terms 0{al) discarded. The various truncation prescriptions 
are extensively studied in ref. [Q. It turns out that the effect on the evolution of e.g. 
F2 is surprisingly large: differences of up to 6% show up at low x ~ 10~^ and are caused 
by terms in NNLO and beyond. The corresponding shift in as{M'^) is estimated to be 
about 0.003. 

• Renormalisation and factorisation scale uncertainties. In Eq. (|I]) and Eq. the renor- 
malisation scale (i?^) and the mass factorisation scale (M^) are both assumed to be equal 
to the momentum transfer Q^. The expressions for the case that and are chosen 
to be unequal can be found in Table gives the shifts in as(M|) when these scales are 
varied independently in the range to AQ^ for different values of a cut made on 
the data. It is seen that the scale dependence is by far the largest contribution to the the- 
oretical uncertainty in and that in particular the mass factorisation scale dependence 
increases strongly with a decreasing cut. 



Summary 



An experimental accuracy of Aq;s(M|) = 0.001-0.002 might be in reach provided the following 
conditions are satisfied: (i) F2 measurements become available in the full HERA kinematic 
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Cut (GeV^) 



gV4 i?' 



:2 = 4g2 



4 

20 
50 



+0.0067 
+0.0032 



-0.0120 
-0.0044 
-0.0029 



-0.0076 
-0.0067 
-0.0061 



+0.0059 
+0.0049 
+0.0042 



Table 3: The theoretical shifts on as{M'^) from scale variations. 



range with systematic and statistical errors of a few percent only; (ii) The dependence of the 
systematic errors on x and is sufficiently well known so that their effects can be absorbed 
in the QCD analysis and (iii) The HERA data can be rehably combined with fixed target F2 
data. 

The various prescriptions of the NLO evolution, which differ by terms of NNLO and beyond, 
cause a theoretical error in ^^(Ml) of about 0.003. If the renormalisation scale {R) and the 
mass factorisation scale (Af) arc varied independently in the range g^/4 to 4g^ a theoretical 
uncertainty on q;s(M|) of about ±0.005 (i?) and ±0.003 (M) is estimated, provided a cut 
of 50 GeV^ is apphed. These errors increase when the cut is lowered. The theoretical un- 
certainties are expected to be reduced substantially once the NNLO splitting functions become 
available. 



I would hke to thank J. Bliimlein, M. Klein and A. Vogt for discussions and comments on 
the manuscript and the DIS97 organizing committee for a stimulating workshop in a beautiful 
location. 
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